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PREFACE 


The work described in this report was performed in 
compliance with the requirements of Contract UI/ARI-88-1, Antenna 
Pattern Study, Task II. Two electromagnetic scattering codes, 
NEC-BSC, Revision 2, and ESP4 were installed on a Compaq 386/20 
computer, together with graphic display codes and auxiliary input 
data codes. The NEC-BSC and ESP4 codes were developed at the 
Electroscience Laboratory of Ohio State University. Several 
revisions were made on these codes by ARI to make them suitable 
for use on the Compaq computer. 
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1.0 INTRODUCTION 


Two electromagnetic scattering codes, NEC-BSC and ESP3, have 
been delivered under a previous contract ( UI/ARI-87-3 ) and have 
been installed on a NASA VAX computer for use by George C. 
Marshall Space Flight Center antenna design personnel. The 
purpose of the work performed under this contract was to update 
the existing codes and certain supplementary software, to install 
the codes on a computer that will be delivered to the customer, 
to provide capability for graphic display of the data to be 
computed by use of the codes and to assist the customer in the 
solution of specific problems that demonstrate the use of the 
codes. With the exception of of one code revision, all of these 
tasks have been performed. 
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2.0 CODE REVISIONS 


Tasks 1 through 3 of the Statement of Work of this contract 
call for acquisition and evaluation of the scattering code NEC- 
BSC, Revision 3. This revision was expected to be released by 
Ohio State University shortly after award of this contract. 
However, it has not yet been released, and recent communications 
with OSU indicate that it will not be released for several more 
months. Therefore, it was not possible to incorporate the 
revised code in this effort. Discussions between all concerned 
parties have led to an agreement that installation of the 
existing version of the code on the PC will be satisfactory. 

No revision of the moment-method code ESP3 was required 
under the contract, inasmuch as the code was not expected to be 
revised. However, use of the code by ARI and by others revealed 
an that an error existed which degraded the quality of the 
computed data under certain conditions. The error was 
subsequently corrected by OSU and a revised version ESP4 was 
released. That version was immediately procured and was 
installed on the Compaq 386/20 computer. 

Several revisions in the two codes were made by ARI and will 
be discussed in detail. 
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3.0 INSTALLATION OF NEC-BSC AND ESP4 ON 386/20 COMPUTER 

The random-access memory requirement for NEC-BSC does not 
exceed that which can be accommodated by the 16-bit Microsoft 
compiler. The requirement is also essentially independent of the 
electrical size (size measured in wavelengths) of the scattering 
body or bodies for which the data are being computed. It depends 
only on the number of scattering elements being modeled, the 
dependency being somewhat linear. it is capable, in its present 
form, of solving scattering problems involving fairly complex 
structures, and is considered adequate for typical NASA uses in 
that respect. 

The ESP4 code, like its ESP3 predecessor, has a very 
different type of memory requirement. The required amount of RAM 
depends on the electrical size of the scattering body. The 
moment method, on which the code is based, involves the formation 
and inversion of a square matrix, the size of which depends on 
the number of "patches" into which the scatterer surface is 
divided. The amount of memory required to perform the required 
computations increases exponentially with model size. An 
assessment of memory requirement for the 386 computer, taking 
into account the type of scattering problems that may be expected 
to be encountered by MSFC design personnel, as well as run time 
and cost, was made by ARI . A compromise was made between these 
factors and the computer RAM was increased to a value of 9 
megabytes. This compromise value is considered adequate to 
accommodate all typical design problems concerning antenna 
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radiating elements and associated electrically conductive 
surfaces. It can be used for scattering bodies which are roughly 
three or four wavelengths (or less) across. 

The large amount of memory required by the ESP4 code cannot 
be accessed by a 16-bit compiler, so that a 32-bit compiler had 
to be acquired. It was also decided that both codes should be 
compiled by use of the same compiler. The initial choice was the 
Lahey compiler. This compiler was procured and used to compile 
the NEC-BSC code, which was then tested on examples foi which the 
correct output data were known. Initial attempts to do the 
compilation were unsuccessful, and it was discovered that the 
compiler would not accept certain statements as formatted in the 
source code. The format of these statements was changed and the 
compilation was accomplished. However, use of the executable 
code on test cases did not produce correct data. Further 
investigation revealed that certain statements in the source 
code, although expressed in good Fortran, could not be properly 
compiled. The vendor was apprised of this deficiency and very 
promptly made a correction to his software, A revised version of 
the compiler was then made available. This version correctly 
compiled the statements mentioned above, and the resulting 
executable code produced correct data when flat-plate scatterer 
models were used. However, when cylindrical scattering models 
were introduced, incorrect data were again generated. 

The errors produced in the executable code were suspected to 
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be related to some optimization feature in the Lahey compiler. 
The source code, as received from OSU, would produce incorrect 
output data when compiled on the VAX in the OPTIMIZE mode. If 
compiled in the NOOPTIMIZE mode, the resulting executable code 
would produce good data. Ohio State personnel experienced the 
same results. It was found by ARI that the source code could be 
correctly compiled on the VAX in the OPTIMIZE mode if a section 
of the main program was removed and made into a subroutine which 
is called at the appropriate point by the main program. This 
change was made, and remains in the cod, as delivered. However, 
compilation on the Compaq PC, using the Lahey compiler, still did 
not produce a good executable code. 

At this point, it was decided that it would not be wise to 
continue the effort to use the Lahey compiler, and that an 
alternate choice should be considered. Acting on that decision, 
the SVS FORTRAN-386 compiler was procured from Science 
International Corporation, Los Altos, California. 

Initial efforts to compile the NEC-BSC source code with the 
SVS compiler were not successful. Again, errors occurred that 
seemed to be associated with optimization. An inquiry was made 
to the vendor and it was learned that the compiler contains five 
optional levels of optimization: 0 through 4. The default level 
is level 4 (the highest), and is the level that was being used. 
A change to level 3 was made and successful compilation was 
accomplished. After being compiled and linked, the resulting 
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executable code was tested and found to generate valid data on 
all test cases on which it was tried. 

After successful compilation and linking of NEC-BSC, the 
ESP4 source code was also successfully compiled and linked, using 
the SVS COMPILER and the Phar Lap linker. It also was found to 
produce valid output data on known test cases. 
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4.0 GRAPHICS CAPABILITY 

A requirement exists for providing a graphic display of the 
data generated by the NEC-BSC and ESP4 codes. These codes 
produce output files which contain data relating to various 
components of the computed field, and for optional pattern 
planes. Thus, the desired graphics code should provide access to 
any polarization component that the user wishes to view and 
should do in convenient manner. It should also perform scaling 
functions which will make most effective use of the plotting area 
and will s lect the most appropriate increments of angle and 
amplitude for the data being displayed. The two graphic codes 
that are provided meet these requirements. 

The first of the two codes, written for use with NEC-BSC, is 
called PLTGTD . The second, for use with ESP4 , is PLTMM. These 
codes make use of a set of graphics routines called GRAFMATICS. 
Both codes request from the user the desired polarization 
component, which he then enters from the keyboard. The code then 
selects the appropriate columns of data from the output file, 
determines whether the pattern cut is of the great-circle or the 
conical type (for NEC-BSC), performs the scaling, sets up the 
grids and displays the data. The resulting plot can then be 
viewed and, if desired, transferred to the printer for hard copy. 
Examples are shown in Figure 1. 
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Figure 1. Example of antenna pattern display. 



5.0 USE OF THE SCATTERING CODES 


The NEC-BSC and ESP4 codes are are stored on the hard disk 
in two subdirectories. NEC-BSC is stored in subdirectory GTD. 
ESP4 is stored in subdirectory MOMENT. They are conveniently 
executed and displayed by use of two batch files, GTD (in the GTD 
subdirectory) and MM (in the MOMENT subdirectory). 

To run the NEC-BSC code, the GTD subdirectory is selected, 
after which the operator simply types the command GTD on the 
keyboard. This command sets up a sequence of actions. The first 
of these actions is to install the GRAPHICS capability. The 
second is to delete any previous plot file. The NEC-BSC 
executable code is then called and executed. In the course of 
execution, it will ask the user the name of the input data file 
to be used in the computation. This file must have already been 
prepared. The output data are then computed and placed in the 
output files. Figure 2 shows the messages that are output to the 
monitor as the code is being executed. One of the output files, 
OUTPUT.DAT, contains the complete output information, including 
input data and other descriptive information. It may be output 
to the printer at any time. The other file, PLOT . DAT , contains 
only numeric data to be accessed by the graphics code PLTGTD . 
This code is then called and requests that the user identify the 
polarization component that he wishes to display. After that 
piece of information has been entered from the keyboard, the 
graphics code selects and displays the requested data. 
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C:\GTD>GTD 


C : \GTD>GRAPHICS 
C:\GTD>DEL PLOT.DAT 
C:\GTD>SCAT 

What is the name of the input file? 
EXAMPLE1.DAT 

ReCE: FAR ZONE PLATE TEST, EXAMPLE 1A. 
ReUN : UNITS IN INCHES 
ReFR : FREQUENCY IN GHZ. 

RePD: PATTERN CUT 
RePG: PLATE GEOMETRY 
ReSG: SOURCE GEOMETRY 
ReXQ: EXECUTE CODE 
ReEN: END CODE 


C:\GTD>PLTGTD 

THE PATTERN IS READY TO BE PLOTTED 


ENTER POLARIZATION COMPONENT DESIRED: 

1 - E-THETA 

2 — E— PHI 

3 - MAJOR AXIS OF POLARIZATION ELLIPSE 

4 - MINOR AXIS OF POLARIZATION ELLIPSE 

5 - TOTAL AMPLITUDE 


Figure 2. Monitor display for execution of NEC-BSC. 
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Execution of ESP4 is accomplished in essentially the same 
way. In this case, the user inputs the command MM from the 
keyboard, rather than GTD. The requests for input file name and 
polarization component appear as before, and the display is 
immediately presented on the monitor. The degree of completion 
is shown on the monitor screen as the impedance matrix is 
computed and subsequently, as the pattern is computed. An 
example of this display is shown in Figure 3. 

Certain changes were made by ARI in the NEC-BSC coce to 
include the E- theta and E-phi components, together with the major 
and minor axes of the polarization ellipse and the total field, 
all in a single plot file. This change permits the plot code 
PLTGTD to access any selected component from the one file. 
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C : \MOMENT>MM 

C : \MOMENT> GRAPH ICS 

C:\MOMENT>DEL PLOT . DAT 

C:\MOMENT>ESP4 
$Enter filename for input 
EXAMPLE1.DAT 
Opening EXAMPLE1.DAT 
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CPU RUN TIME FOR RUN 1 GEOMETRY 1 = 70.25 SECONDS 


TOTAL CPU RUN TIME * 70.53 SECONDS 


Programmed STOP 
C:\MOMENT>PLTMM 

THE PATTERN IS READY TO BE PLOTTED 


ENTER POLARIZATION COMPONENT DESIRED: 

1 - E-THETA 

2 - E-PHI 


Figure 3. Monitor Display for ESP4. 
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6 . 0 SPECIFIC STUDIES 


6.1 SCATTERING STUDY FOR SHUTTLE EXTERNAL TANK 

During the course of this contract, MSFC antenna design 
engineers expressed a concern with the type of scattering that 
will occur from the external fuel tank of the Shuttle when radar 
transponder antennas are located on the forward skirts of the 
solid-rocket boosters (SRB). The NEC-BSC code was viewed as a 
useful tool in predicting the composite radiation pattern that 
would result from the combination of direct radiation and the 
radiation that would be scattered by the external tank. However, 
a question existed regarding the validity of the data that would 
be derived by use of the code. The external tank would be 
modeled as a smooth, perfectly-conducting right circular 
cylinder. That model does not exactly agree with the actual 
shape of that part of the tank that lies alongside the antenna 
location. The surface of the tank in that region (the inter-tank 
section) does not have a smooth cylindrical surface. Additional 
strength is provided in the intertank structure by the use of 
external longitudinal stringers. These stringers produce a 
somewhat corrugated shape at the scattering surface of the tank. 

A study was performed by ARI to assess the validity of the 
data that would be derived by use of the NEC-BSC code, treating 
the intertank skin as a smooth cylindrical surface. The work 
performed in that study and the conclusions that were drawn are 
described in Report ARI/R-89-573 , ELECTROMAGNETIC SCATTERING BY 
SHUTTLE EXTERNAL TANK SURFACE, dated February 27, 1989. 
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6.2 OBSCURATION COMPUTATION FOR SRB RADAR ANTENNAS 

An intermediate step in determining the adequacy of the 
radar transponder antenna patterns for two specified radar 
stations was computation of the optical shadow produced by the 
various bodies comprising the shuttle cluster in the launch 
configuration for a source located at an antenna location. This 
computation was performed by ARI , using the computer code SHADOW. 
In the absence of scale-model antenna pattern measurements, the 
results of the computation were used by NASA in preliminary 
assessment of tne adequacy of the antenna arrangement and the 
radar station layout. A detailed description of the obscuration 
computation is given in report ARI/R-89-575, TRANSMISSION PATH 
OBSCURATION FOR RADAR TRANSPONDER ANTENNAS ON SHUTTLE SRB, dated 
March 24, 1989. 
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7.0 CONCLUSIONS 


The two scattering codes NEC-BSC and ESP4 have been 
installed on a COMPAQ DESKPRO 386/20 computer, together with 
graphics software and auxiliary input data codes. All required 
revisions have been made and the codes are ready for use. It is 
believed that the codes, in their present form, will provide a 
valuable predictive capability for antenna patterns associated 
with large and complex scattering bodies with which NASA is 
concerned. The absence of the expected second revision of NEC- 
BSC code is regrettable, but the current revision, installed on 
the Compaq computer, with graphic display, is expected to be 
highly useful. It is recommended that the third revision be 
acquired and installed on the computer when it becomes available. 
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APPENDIX A 
INPUT DATA CODES 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


PROGRAM BSCINPUT 

THIS PROGRAM PROVIDES AN INTERACTIVE MEANS FOR ENTERING INPUT 
DATA INTO THE OHIO STATE UNIVERSITY NEC-BSC PROGRAM. THE 
REQUESTED PARAMETER VALUES ARE WRITTEN TO A DATA FILE 
TO BE CALLED BY PROGRAM "SCAT". 

THIS PROGRAM WAS WRITTEN BY: 


J. WARREN HARPER 
APPLIED RESEARCH, INC. 
5025 BRADFORD BLVD. 
HUNTSVILLE, ALABAMA 35805 


Change section written by Mike Vinson 


Revision date: January 27, 1989 


Character *3 
Character *36 
Character *36 
Character*2 
Character *1 
Character *30 
Character*40 
Character *40 


STRING 

Temp 

INPUT 

NZERO 

Answer 

Descr iption( 33 ) 
Line 

Section ( 10) 


Description(l) 
Description 2) 
Description ( 3 ) 
Description( 4 ) 
Description ( 5 ) 
Description^) 
Description ( 7 ) 
Description ( 8 ) 
Description ( 9 ) 
Description ( 10 ) 
Description( 11 ) 
Description (12) 
Description (13) 
Description ( 14 ) 
Description ( 15 ) 
Description (16 ) 
Description (17 ) 
Description 18 ) 
Description ( 19 ) 
Description 20) 
Description 21 ) 
Description(22) 
Description(23) 
Description ( 24 ) 
Description(25) 
Description! 26 ) 
Description(27) 
Description(28) 
Description 29 ) 
Description 30 ) 
Description ( 31 ) 
Description ( 32 ) 
Description ( 33 ) 


'BP: Back or Bistatic Scatter' 
'CE: Last or only Comment' 

'CG: Cylinder Geometry' 

'CM: Comment Card' 

'EN: End program' 

'FM: Swept Frequencies' 

'FR: Frequency' 

'GP: Infinite Ground Plane' 

'LP: Line Printer Output' 

'NC: Next set of Cylinders' 

'NG: No Ground Plane' 

•NP: Next Set of Plates' 

'NR: Next Set of Receivers' 

'NS: Next Set of Sources' 

'NX: Next Problem' 

'PD: Far Zone Pattern Cut' 

'PG: Plate Geometry' 

'PN: Near Zone Pattern Cut' 

•PP: Plotter output' 

•PR: Gain or Coupling Factors' 
•RA: Receiver Array Geometry' 
'RD: Far Zone Range' 

•RG: Receiver Geometry' 

'RM: NEC - MM Receiver Input' 
'RT: Rotate-Translate Geometry' 
'SA: Source Array Geometry' 

'SG: Source Geometry' 

'SM: NEC - MM Source Input' 

•TO: Test Options' 

'UF: Model Scale Factor' 

'UN: Units of Geometry' 

’US: Units of Source Size' 

'XQ: Execute Code' 
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c 

Write 0,300) 

300 Format ( 1H$, ’ Would you like to edit the existing file? ’ ) 

Read <*,310.> Answer 

310 Format(Al) 

320 Format (A40) 
c 

c See if the user wants to edit the current file 

c 

I f »; answer .eq. ’ V ' .or. answer .eq. 'y') then 
c 

c Enter the change section 

c 

OpenC Unit = 1, File = ' Input. bak' , Status = ’Old’) 

Open ( Unit = 2, File = dnput.dat’, Status = 'Unknown') 

Read(l, 320, End =999) Line 
400 Call Sear ch (Line, Description, k > 

Sec t ion d) = Line 
n = 0 

i count = 1 

c Do While (n .eq. 0) 

405 Continue 

Read < 1 , 320, End =999) Line 
Call Search (Line, Description, n) 
i count = i count + 1 
Sec t ion d count ) = Line 

If< n .eq. 0) goto 405 
c End do 

i c o u n t = i c o u n t - 1 

Cal 1 Change (Sec tion, k, Description, i count) 
goto 400 

999 Continue 

Cal 1 ChangetSection, k, Description, i count) 

C 1 ose ( 2 > 

Else 

□ PEN d, FI LE = ' F0R005.DAT' , STATUS = ’ NEW’ ) 

C 

STRING = ’CE:’ 

WRITE C*, 1 > 

1 FORMATdH*, ’ Enter Comment (max 36 chars. 1 ’ > 

READO, 2) INPUT 

2 FORMAT (ASS) 

WRITE (1,3) STRING, INPUT 

3 FORMAT (A3, ASS) 

Get units for this run from the user. 

1000 WRITE (*, 200) 

200 FORMATdX, F Enter number indicating units to be used -'/15X, 

.'1 = Meters' /15X, f 2 = Feet'/15X,'3 = Inches'/) 

READ(*, 2) INPUT 

If( Input .gt. '3' .or. Input .It. '1' ) Then 

Write(*,*) 'Your selection was out of range. Please try again 
goto 1000 
End if 
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c 

STRING*’ UN: » 

WRITE(1,9) STRING 
WRITE (1,2) INPUT 
C 

C Get frequency from user. 

C 

WRITEC*, 5) 

5 FORMATC IX, 'FREQUENCY IN GHZ.?’) 

READ*:*, 2) INPUT 
STRING*’ FR: ' 

WRITE*. 1,9) STRING 
WRITE Cl, 2) INPUT 

7 FORMAT*: 12) 

8 FORMAT C I 1 , 1X,A1) 

9 FORMAT (A3) 

C 

WRITE*:*, 10) 

10 F0RMATC1X,' *** Orientation of Pattern Axes ***'/) 

STRING*' PD: ' 

WRITEC*, 15) 

15 FORMAT (IX, 'THETA, PHI FOR Z AXIS, THETA, PHI FOR X AXIS (1 LINE* 
RE ADC*, 2) INPUT 
WRITE*: 1,9) STRING 
WRITE (1,2) INPUT 
C 

c *** TYPE OF PATTERN CUT *** 

C 

WRITE*:*, 210) 

210 FORMAT (IX, 'TYPE OF PATTERN CUT DESIRED (GREAT CIRCLE OR CONICAL)’ 
1/' T = CONICAL CUT (CONSTANT THETA)'/' F = GREr 

2 CIRCLE CUT (CONSTANT PHI)’//’ ENTER T OR F, FOLLOWED BY VALUE O 1 
3THE FIXED ANGLE, ON ONE LINE.') 

READ*:*, 2) INPUT 
WRITEC 1,2) INPUT 
C 

C *** ANGULAR RANGE DESIRED **★ 

C 

WRITE*:*, 220) 

220 FORMAT (IX, 'ANGULAR RANGE DESIRED FOR PATTERN: INITIAL ANGLE, FIN f 

1 ANGLE, ANGULAR INCREMENT’!) 

READ*:*, 2) INPUT 
WRITEC 1,2) INPUT 
C 

C *** PLATE GEOMETRY *** 

C 

WRITE*:*, 40) 

40 FORMAT (IX' HOW MANY PLATES?’) 

READC*, *) NOPLT 
IFCNOPLT.LT. 1) GOTO 82 
C 

STRING*’ PG: ' 

DO 80 MP*1, NOPLT 
WRITE (1,9) STRING 
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C 

WRITE<*,50) MP 

50 FORMATdX, ' PLATE NUMBER ',12,/' HOW MANY CORNERS?' > 

READ(*, 7) NUM 
NZER0='0» 

WRITE< 1,8) NUM, NZERO 
C 

DO 70 ME =1, NUM 
WRITEC*, 60) ME 

SO FORMAT (IX, ' X,Y, Z POSITION OF CORNER NUMBER ',12) 

READ (*, 2) INPUT 
WRITEC 1,2) INPUT 
70 CONTINUE 
80 CONTINUE 

*** CYLINDER GEOMETRY *** 

82 WRITE*:*, 84) 

84 FORMAT*: IX, 'HOW MANY CYLINDERS?') 

READ*.*,*) NCYL 
IFCNCYL. LT. 1) GOTO 99 
DO 94 N=1 , NCYL 
STRING='CG: ' 

WRITE*: 1,9) STRING 
WRITE*:*, 8S) N 

86 FORMAT*: IX, ' LOCATION *.X,Y,Z) OF ORIGIN, CYLINDER NO. ',11) 
READ(*,2) INPUT 
WRITE*: 1,2) INPUT 

WRITE*.*, 88) N 

88 FORMAT*: IX, 'CYLINDER NO. ',11,' OR I ENTAT I ON' / 1 X , ' THETA -Z , PHI-Z, 
1ETA-X, PH I -X ' ) 

READC*, 2) INPUT 
WRITEC 1,2) INPUT 

WRITE*:*, 90) 

90 FORMAT*: IX, ’ CYLINDER RADII: RX, RY') 

READC*, 2) INPUT 
WRITE (1,2) INPUT 

WRITE*:*, 92) 

92 FORMAT ( 1 X, ' POSITIONS S< ANGLES OF END CAPS' / 1 X ,’ POSI TION & ANGLE 
1R NEGATIVE END, POSITION & ANGLE FOR POSITIVE END.') 

READ*:*, 2) INPUT 
WRITE (1,2) INPUT 
94 CONTINUE 


*** SOURCE GEOMETRY *** 

99 WRITE (*, 100) 

100 FORMAT < IX,' *** SOURCE GEOMETRY ***'/) 

WRITE(*, 110) 

110 FORMATdX, 'HOW MANY SOURCE ELEMENTS?') 

READ ( * , * ) NWIRES 
STRING 3 ' SG: ' 
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c 

DO 180 MS=1 , NWIRES 
WRITE (1,9) STRING 
WRITE (*, 120) MS 

120 FORMAT (IX, ’ SOURCE ELEMENT NUMBER' , 12 /) 

WRITE (*, 130) 

130 FORMAT (IX, ' X , Y, Z POSITION OF ELEMENT CENTER?') 

READ<*,2) INPUT 
WRITEC 1, 2) INPUT 
C 

WRITE*:*, 140) 

140 FORMAT Cl X, 'THETA PHI ANGLES FOR LENGTH !< WIDTH VECTORS'/' 

1 C THETA -L, PHI-L, THETA-W, PHI-W') 

READC*, 2) INPUT 
WRITEC 1,2) INPUT 
C 

WRITE*:*, ISO) MS 

ISO FORMAT*: IX, ' TYPE, LENGTH AND WIDTH OF ELEMENT' , 12/ ' TYPE:', 

1 -1 UNIFORM CURRENT DISTRIBUTION'/' -2 

2PIECEWISE SINUSOIDAL DISTRIBUTION'/) 

READ*:*, 2) INPUT 
WR I TE (1,2) INPUT 
C 

WRITE*:*, 170) 3 

170 FORMAT C 1 X , ' EXl ITATION (MAGNITUDE, PHASE) FOR ELEMENT ', 1 2 > 

READ <*, 2) INPUT 
WRITE (1,2) INPUT 
180 CONTINUE 

STRING = ' XG: ' 

WRITE Cl, 9) STRING 
STRING=’ EN: ' 

WRITEC 1,9) STRING 
End if 
CLOSE Cl) 

END 


Subroutine Search (String, Array, n) 
c 

Character*BO String 
Characterise Array (33) 
c 

c Search for the first two letters of the string in the first two 

c letters of the elements of the array, 

c 

Do 10 i = 1,33 

IfC string*: 1:2) .eq. Ar ray C i ) C 1 : 2) ) then 

n = i 
Return 
End if 

10 Continue 
n = 0 
Return 
End 
c 
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c 


Subroutine Change (Sect ion, k, Description, i count) 


c 


c 


c 


c 

c 


10 


90 


c 


300 


80 


1 00 


c 

c 

c 

c 

c 


500 


Character*40 

Character*40 

Character*30 

Characters 

I nteger*4 
I nteger*4 
I nteger*4 
I nteger *4 

Character*l 
I nteger*2 


Section( 10) 

Line 

Descript ion (33) 
Answer 

k 

1 count 
Offset 
jof f set 

Btemp<2) 

I temp 


Equivalence ( I temp, Btemp(l) ) 


Write(*, *) 

Write(*, *> 

Write(*,*> 'Command card found: ' , Description(k) < 1 : 3) 

Write<*,*> 'This is described as: ' , Description(k) (4: 30) 

Write<*, *) 

Write!*,*) 'These are the cards in this set:' 

Do 10 i = 1,1 count 

Write*;*, *) Section(i) 

Continue 
Write*.*, *) 

Write*.*, 90) 

Format *.' 1H$, ’ Do you want to change any of the cards in this SET 

Read (*, 1 10) Answer 

Write*;*,*) 

Write (*,*) 

If( Answer .eq. 'Y' .or. Answer .eq. 'y') Then 
Do 20 i = 1,1 count 
Write*;*,*) Section(i) 

Write*:*, SO) 

Format *. / 1 X ,' Note: To insert a new card',/ 

' before this card, just',/ 

' type INSERT To delete',/ 

' a card, type DELETE.'/) 

Write*:*, 100) 

Format < 1H$, ' Do you want to Change this card? ' > 

Read*:*, 120) Line 
Answer = Lined: 1) 


To test all possible combinations of the word insert in both 
upper and lower case, convert each of the first six letters 
to upper case. 

Do 500 k = 1,6 

Btemp(l) = Line(k:k) 

if( itemp .ge. 97 .and. itemp . le. 122 ) then 
I temp = Itemp - 32 
Line(k:k) = Btemp(l) 

End if 
Continue 
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70 


125 


130 

200 

20 

30 

110 

120 


If( Line(l:6) .eq. ’INSERT’ ) then 
Call Insert 
Goto 300 
End if 

iff. Lined :6) .eq. ’DELETE' ) Goto 20 
If< Answer .eq. ’Y’ .or. Answer .eq. ’y’) Then 
Write(*, * ) 

Write(*,70> 

Format (// 1 X , ’ Note: Type to begin and end inserted dat 

’ for this card. Type # to delete' , / 

’ characters.’/ ) 

Write*!*, *> 'Enter the New card below (Space = no change 
Write (*,*) Sec tion C i ) 

Read (*, 120) Line 
Ins_flag = 0 
Offset = 0 
Do 200 j = 1 , 40 

I f ( Line(j:j) .eq. ’ ’ ) then 

Ins_flag = 1 - Ins_flag 
Offset = Offset + 1 
Goto 200 
End if 

If( Line(j:j) .eq. ’#’) then 
Do 125 m = j -Off set, 39 

section< i ) <m: m!> = sec tion ( i ) Cm+1 : m+1 ) 

Continue 

Offset = Offset + 1 
Goto 200 
End i f 

If( Line(j:j> . ne. ’ ’ .and. lns_flag .eq. 0) then 

joffset = j - Offset 

Sec tion ( i )( joff set : jof f set > ~ line(j:j) 

Else if( ins_flag .eq. 1 ) then 
Do 130 m = 39, j-Of fset, -1 

Sectiont i ) (m+1 : m+1 > = Sec t ion ( i ) (m: m) 

Continue 

joffset * j - Offset 

Section(i) (joffset: joffset) = line(j:j> 

End if 
Continue 
End if 
Write*!*,*) 

write (2, 120) Section(i) 

Continue 

Else 

Do 30 i = l,i count 

write (2, 120) Section(i) 

Continue 
End if 
Format (A1 ) 

Format ( A40) 

Return 

End 
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<N CO 


c 

c 


Subroutine Insert 


Character*40 Line 
Character*l Answer 

c 

1 Write(*,10) 

10 FormatC/’ Type the card to be inserted:’) 

Read (*, 100) Line 
Wnte(*,200) Line 
WriteC*, 300) 

Read (*,90) Answer 

I f C Answer .eq. 'Y’ .or. answer .eq. ’y’) then 
w r i t e ( 2 , 1 00 ) Line 
Else 

Goto 1 
End if 

90 Format (Al) 

1 00 For ma t ( A40 ) 

00 Format (//IX, A40/ ) 

00 Format ( 1H$, ’ Is this card correct? ’) 

Return 

End 
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PROGRAM INDAT4 


C This program provides an interactive means for entering input 

C data into the ohio state university esp4 program. The requested 

C parameter values are written to a data file* specified by the 

C user, in the proper order to be called by esp4. 

C 

C This program was written by: 

r 

C J. WARREN HARPER 

C APPLIED RESEARCH, INC. 

C 5025 BRADFORD AVE. 

C HUNTSVILLE, ALABAMA 35S05 


c 

c Modified by: Mike Vinson 

c 


C 


C 


c 


c 


DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 

DIMENSION 


NCNRS( 10) 

X!10). Y ( 10 ) , Z(10) 

NPLA( 10) , BDSK! 10), VGA< 10), ZLDA( 10) 

P C N < 3 , 10, 10), IA( 10), IB ( 10), SEGM ( 10) 
IREC(IO), IPN(IO), IGS(IO), NAS (10) 

NSA < 10) , VG( 10) , ZLD( 10) , IFMC 10) 
IABFP(IO), IABAP(IO) 

VLG( 10), ZL( 10) 


CHARACTER *1 USEWRS, GEOPRT, DFLT , FARZN, USEPLT, MCOUPL, WRTIMP 
CHARACTER *1 CHFREQ, CHOP AT, CHGRAD, CHGCON, Chang e 
Charac ter *40 Filename, Variable 
Character*! Answer 


COMPLEX VLG 
COMPLEX ZL 
COMPLEX VGA 
COMPLEX ZLDA 
Complex*8 Zsht<30) 


10 

20 

30 


WRITE ! *, 10) 
WRITE!*, 20) 

WR I TE ( * , 30 > 
FORMAT!' This 
FORMAT! ' 
FORMAT ( ' ' ) 


program !ESP4) is based on the method 
*********************** ' ) 


of moments' 1 ) 


WRITE!*, 2000) 

2000 FORMAT! 1H$, 'Do you wish to change an existing input file? (Y/N) 


READ!*, 50) CHGhIL 

IF! CHGFIL . EG. 'N' .OR. CHGFIL . EQ. 'n' ) GOTO 8 
Write!*, 201 1 ) 

2011 Format ( li-tt, 'Enter the filename to edit: ') 


Read!*, 431) Filename 
431 format (A40) 

Open! Unit=iO, F i 1 e=Fi 1 ename > 3tatus= 'old ! ) 
ICHG=0 

£ 

C READ IN VALUES FROM EXISTING DATA FILE 
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READ(10, *) NGO, NPRINT, NRUNS, NWGS, IWR, IWRZT, INT, INTP, 
I NWR , IRGM, IFIL, RF, INDZI 
If( Indzi. . n e . 0) Then 

Read <10, # ) Fmcl, Fmc2, Dfzi, Dff, Irsi2, Thrd, Phrd, 
Thine, Rhine 

End if 


I f ( Indzi . eq. 

0) 

Th en 



READ ( 10,*) 

IFE, 

I PFE, 

FNDFE, 

PHFE 

READ( 10, *> 

IFA, 

IPFA, 

FNDFA, 

THFA 

READ( 10, *) 

ISE, 

IPSE, 

FNDSE, 

PHSE, THIN, PHIN 

READ ( 10, *> 
End if 

ISA, 

IPSA, 

FNDSA, 

THSA 


INTD, 


READ (10, *> FMC, CMM, A 

c IF (USEPLT . ne. ' N ' .OR. USEPLT . ne. 'n'> Then 

READ (10, * ) NPLTS 
If( Nplts . ne. 0) Then 
DO 1520 NPL* 1, NPLTS 

READ (10, * ) NCNRS(NPL), SEGM(NPL), IREC(NPL), IPN(NPL), 
IGS(NPL), Zsht(Npl) 

DO 1510 NCMR=1 , NCNRS(NPL) 

READ (10, *) PCNd, NCNR, NPL). PCN ( 2, NCNR, NPL ) , 

PCN(3, NCNR, IMPL) 

1510 CONTINUE 

1520 CONTINUE 

End if 

c 1530 IF ( USEWRS. EQ. 'N OR. USEWRS. F.Q. 'n'l GOTO 310 
READ (10, * > IWRZM, IRDZM 

If ( Jnwr . ne. 0 . and. Irgin . ne. 0) Then 
READ( 10, *> NM, NP, MAT, NFPT , NFS1, NFS2 

C 

DO 1550 1=1, NP 

’ READ( 10, *> X( I ), Y( I ), Z( I ) 

1550 CONTINUE 

C 

DO 1560 1=1, NM 

READ ( 10, #) I A < I ) , I B ( I ) 
uir i be Of-, ■* ) ia(i), ib(i) 

1560 CONTINUE 

C 

DO 1570 1 = 1, NFPT 

READ ( 10, *> IFM ( I ) , I ABFP ( I ) , VLG ( I ) , ZL(I) 
urite<*, »•> IFM ( I ) , IABFP(I), VLG(I), ZL(I) 

1570 CONTINUE 

C 

DO 1580 1=1, NAT 

READ ( 10, *) NAS ( I ) , IABAP ( I ) , NPL A ( I ) , VGA ( I ) , ZLDA ( I ) , BDSK ( I ) 
wr ite(*, *> NAS ( I >, IABAP (I ), NPLA( I ) , VGA ( I ) , ZLDA ( I ) , BDSK ( I ) 

1580 CONTINUE 

End i f 
REWIND 10 
r 

C END OF READING OF DATA FILE 


Variable = 'Frequency' 

Call Rchange ( Variab 1 e, FMC, 0) 

3200 FORMAT ( / ) 

C 

c 2002 IF ( USEWRS. EQ. 'N OR. USEWRS. EQ. 'n') GOTO 2006 

26 


I f < Nwg s . g t . 0 ) I hen 

Variable = 'Wire Radius' 

Call Rchange( Variable. A. 0) 

r 

Variable = 'Wire Conductivity' 

Call Rchange( Variable. Cmm, 0) 

End if 
c 

If! INDZI . ne. 0) then 

5020 Wr i t e (#. # ) 'The current values for Frequency Sweep are: ' 

wr i te ( #, 5010 ) Fmcl, Fmc2» Dfzi. Dff. Irsl2. Thrd. Phrd. 

Thri, Phri 

5010 Format!5x> '1 -Begin of sweep ( MHz ) ' . FS. 2. / 

5x. '2 —End of sweep ( Mh z ) ' . F8. 2. / 

5x» '3 -Step size for Imdepance Matrix ( Mh z ) ' , F8. 2. / 

5x. '4 -Step size for MM computation ( Mh z ) ' . F8. 2. / 

5x. '5 -Indicator of radiation or scattering pr ob 1 em ' , 12, / 
10x» ' 1 - Implies radiation problem',/ 

lOx, ' 2 - Implies scattering problem',/ 

lOx, '-2 - Implies scattering including incident wave',/ 

5x, '6 -Angle for radiated or scattered Pei Id (Degrees)',/ 

1 0 x , 'Theta = ',F8. 2, 5x, 'Phi = ',F8. 2, / 

5x,'7 -Angle of incident for scattering problem (Degrees)',/ 
lOx, 'Theta = ' , FB. 2, 5x , 'Ph i = ',F8. 2//> 

Wr i t e < 6050 ) 

Read (a-, 50) Answer 

If (Answer . eq. 'y* . or. Answer . eq. 'V') Then 
Write!#, 2035) 

Read <#,«-) I 
I f < I eq. 1 ) then 

Variable = 'Beginning of Sweep (MHz)' 

Call Rchange! Variable, Fmcl, 0) 

Else if < I . eq. 2) then 

Var iab 1 e = ' End o f Sweep ( MH z ) ' 

Call R change! Variable, Fmc2, 0) 

Else if(I eq. 3) then 

Variable = 'Step Size for Z matrix (MHz)' 

Call R change! Variable, Dfzi, 0) 

Else if( I . eq. 4) then 

Variable = 'Step Size for MM Computation (Mhz)' 

Call Rchange! Variable, DFF, 0) 

Else if! I . eq. 5) then 

5025 Variable = 'Radiation or Scattering Problem' 

Call Ichange (' Variab le, IRS12, 0) 
j = I r s 1 2 

if! j. ne. 1 . or. j. ne. 2 . or. j. ne. -2 ) then 

Write!#,#) 'Input out of range. Try again' 

Goto 5025 
En d if 

Else if! l.eq. 6) then 
Variable = 'Theta' 

Call Rchange! Variable, Thrd, 0) 

Variable = 'Phi ' 

Call Rchange! Variable, Phrd, 0) 

Else if! I.eq. 7) then 
Variable = 'Theta' 

Call Rchange! Variable, Thri, 0) 

Var i ab 1 e = 'Phi' 
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Call R channel variable* Fhri, 0) 

End if 
Goto 5020 
End if 
Else 

Wr i te ( *, 6020 ) If e, Ipfe, Fndfe, Phfe 
6020 Format( ' The current values of Pattern Data are: ',// 

5x> '1 -Compute near or far zone elevation radiation pattern' 

, 12/, 5x, '2 -Write elevational pattern to plot file'. 12* / 

5x« '3 -Elevation angle increment (Deg) '., FS. 2, / 

5x* '4 -Constant phi elevational angle (Deg ) ' » F8. 2/ ) 
c 

Write!*, 6030) IPs. Ipfa, Fndfa, Thfa 
6030 Format(5x> '5 -Compute near or far zone azimuth radiation pattern' 
* 12/* 5x* '6 -Write azimuthal pattern to plot file'* 12/# 

5x* '7 -Azimuth angle increment (Deg) ',F8. 2/. 

5x# '8 -Constant theta azimuthal angle ( Deg ) ' , F8. 2/ ) 
c 
c 

Write<*, 6050) 

6050 Format(lH$, 'Would you like to change any of these? ') 
Read(*»50) Answer 


c 


6060 

6070 


c 


6080 


c 


6040 


Wr i te ( *, 6060 ) 

Format!//' Scattering Section: '/) 

Wr i te ( *, 6070 ) Ise, Ipse# Fndse* Phse 

Fdrmat(5x, '1 -Type of Scattering in Elevation Plane: ';I2#/10x» 


5 x 
5 x 
5 x 


' 1 
'2 
'3 


-Do not compute backscatter' 

-Compute Backscatter pattern 
-Compute Bistatic scattering 
-Computer Forward scattering 
-Write scattering pattern to 
-Angle increment for scattering ';F8. 2#/ 
-Constant Phi angle (Deg) for scattering 


/ lOx < 

* / lOx ; 

pattern ', /10x# 
pattern ' * / 

Plot f ile 12/, 


F8. 2/) 


Wr i t e ( *, 6080 ) Isa, Ipsa, Fndsa, Thsa 

Format(5x, '5 -Type of Scattering in Azimuthal Plane: ',12,/ 

5x, '6 -Write scattering pattern to Plot file', 12/, 

5x, '7 -Angle increment for scattering ' , F8. 2, / 

5x, '8 -Constant Theta angle (Deg) for scattering '» F8. 2/ ) 


Wri te < *, 6040 > Thin, Phin 

Format(5x, '9 -Theta angle (deg) of incident wave in Bistatic', 
' Scatter ing ', F8. 2, /, 

4 x , '10 -Phi angle (deg) of incident wave in Bistatic'# 

' Scatter ing ', F8. 2// ) 

Write!*, 6050) 

Read!*, 50) Answer 
End if 


2009 WRITE!*, 2010) 

2010 FORMAT < 1H$, 'Would 
READ(*, 50) Answer 
IF ( Answer . eq. ' V ' 

Wr i t e ( *, * > ' 

Wr i te ( *, * ) ' 

Write!*,*) ' 
wr i te < *, * ) 

Wr i te (*> 2720) 
2720 Format ( 1H$, 'Your 


you like to change the Plate Data? (Y/N) ') 

'y') Then 


OR. Answer . eq. 

Option Menu ' 

Add a plate' 
Change a plate' 


1 - 
n _ 


choice: 
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Read < *, * ) I 
If < I . eq. 1 ) Then 
NPLTS = NPLTS + 1 
NPL = NPLTS 
Ichg = 1 

CALL PLATE ( NPL, NCNRS, SEGM, IREC, IPN, IGS, PCN, ICHG, 
Zsht) 

Else 

Write ( *> 2030 ) Nplts 

2030 Forma t ( 1 H$ ' Th ere are now', 12, ' plates. ') 

Write!*, 2035) 

2035 Format ( 1H$ 'Wh i ch would you like to change? ') 

READ ! *, * ) NPL 
Ichg = 0 

CALL PI. ATE ( NPL, NCNRS, SEGM, IREC, IPN, IGS. PCN, ICHG, 
Zsht) 

End if 
End if 


2005 WRITE.!*. 2060) 

2060 FORMAT ! 1H$, 'Would you like to change the Wire Data? (Y/N) ' ) 

READ < *, 50 ) Answer 

IF< Answer eq. 'Y' OR. Answer . eq. 'y') l hen 
I c h q = 0 

CALL WIRE < NM» NP, NAT, NFPT, NSF1, NSF2, ICHG) 

End i f 
GOTO 010 

****•**•»•** END OF CHANGE SECTION; BEGIN ORIGINAL INPUT ********* 


B I CHG~ 1 

Wr i te ! *> 310) 

310 Format! 1H$, 'Enter the filename to write: ') 

Read (*. 431 > Filename 

Op en ( Un i t- 1 0, F i 1 e=F i 1 ename » Sta t u 5= ' new ' ) 

WRITE!*, 12) 

12 F0RMAT!1H$, 'Do you wish to specify any of the following', 

' parameters? (Y/N) '/) 

WRITE ( *, 14) 

14 FORMAT! 10X, '* Number of runs to be made'/lOX, 

'* Print of modal currents '/10X, 

'* Print of impedance matrix'/lOX, 

'* Define frequency sweep '/10X, 

'* Number of Simp sons-ru 1 e intervals for: '/15X, 
' Wir e-to— wir e imp edanc es ' / 1 5X» 
'Surface-patch monop o 1 es '/ 15X, 

'Disk monopoles'//) 


READ < *, 50 ) DFLT 

IF(DFLT. EQ. ' Y ' . OR. DFLT. EQ. 'y') Then 
Variable = 'Number of Runs' 

Call Ichanget Variable, Nruns. Ichg) 

Variable = 'Number of Wire Geometries per Run' 
Call Ichanget Variable, Nwgs, Ichg) 
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a cj o o o f o o 


WRITE ( *, i 025 ) 

1025 FORMAT (/IX, MODAL CURRENT PRINTOUT: '/6X, 

' 0 NO MODAL CURRENT PR INTOUT ' /6X, 

'1 MODAL CURRENTS PLUS WIRE/PLATE GEOMETRY PRINTED ' ) 
Variable = 'Modal Current Printout' 

Call Ichange( Variable« Iwr, Ichg) 

r 

WR I TE ( *, 1030) 

1030 FORMAT (/IX/ 'WRITE IMPEDANCE MATRIX IN OUTPUT FILE? (Y/N)') 

READ ( *, 50 ) WRTIMP 

C 

IF ( WRTIMP. EG. 'Y'. OR. WRTIMP. EG. 'y') Then 
1032 IWRZT=1 

Else 

I WRZT=0 
End if 
C 

WRITE (r,, 1040) 

1040 FORMAT ( / 1 X > 'NUMBER OF SIMPSONS-RULE INTEGRATION INTERVALS: '// 

IX. 'FOR WIRE-TO-WIRE IMPEDANCES?') 

READ(*. *) I NT 
WRITE ( *, 1050) 

1050 F0RMAT( /I X, 'FOR SURFACE-PATCH MODULES?') 

READ!*,*) I NTP 
WR I TE ( *> 1060) 

1060 FORMAT (/IX. 'FOR DISK MONOPOLES?') 

READ<*. *) INTO 
WRITE(*. 1070) 

1070 FORMAT (//IX, 'HOW IS THE WIRE GEOMETRY DEFINED? '/ 10X, 

'0 BY SUBROUTINE WGEOM'/lOX, 

'1 DEFINED VIA THE INPUT FILE') 

READ ( *, * ) IRGM 
Else 


APPLY DEFAULT VALUES 

READ 1: 

NRUNS*i 

NWGS=1 
I WR=0 
IWRZT=i 
INT=4 
INTP=6 
INTD=1S 
I RGM= 1 
End if 


INPUT NON-DEFAULT INFORMATION 

35 WRITE (*,40) 

40 FORMAT (/IX, 'DOES THE MODEL CONTAIN WIRES? (Y/N)') 
READ (*,50) USEWRS 
WRITE!*, 45) 

45 FORMAT (/IX, 'DOES THE MODEL CONTAIN PLATES? (Y/N)') 
READ ( *, 50) USEPLT 

50 FORMAT (Al) 

IF ( USEWRS. EQ. ' Y '. OR. USEWRS. EQ. 'y') INWR = 1 


v| M 


WRITE ( *, 60) 

60 FORMAT (/IX, 'SPECIFY REQUIRED OUTPUT: '/10X. 

'1 PRINT WIRE & PLATE GEOMETRY '/ 10X, 

'2 PRINT INPUT PARAMETERS & WIRE/PLATE GEOMETRY '/ 10X 
'3 PRINT NOTHING') 

C > 

READ**,*) NPRINT 

C 

WRITE**, 70) 

70 FORMAT (/IX, 'DO YOU WANT TO PRINT INPUT DATA BEFORE COMPUTING? ' ) 
READ**, 50) I FIRST 

IF* IFIRST. EQ. 'Y '. OR. IFIRST. EQ. 'y') THEN 
NG0=0 
ELSE 
NG0=1 
END IF 

Q 

Write**, 710) 

710 Format * 1H$, 'Do you want a Frequency Sweep? ') 

Read**, 50) Answer 

If* Answer . eq. 'Y' . or. Answer eq. 'y') then 
Write (*, 720) 

720 Format* ' i -Frequency sweep with impedance matrix ', 

' interpolatibn method',/' 2 -Frequency sweep with ', 
'"improved" impedance matrix in ter p o la t i on method') 
Write**, 725) 

725 Format* 1H$, 'Which would you prefer? ') 

Read**,*) INDZI 
c 

Variable = 'Beginning of Sweep (MHz)' 

Call Rchange* Variable, Fmcl, 1) 

Variable = End of Sweep (MHz)' 

Call Rchange* Variable, Fmc2, 1) 

Variable = 'Step Size for Z matrix (MHz)' 

Call Rchange* Variable, Dfzi, 1) 

Variable = 'Step Size for MM Computation (Mhz)' 

Call Rchange* Variable, DFF, 1) 

30 Write**, 735) 

35 Format*' Indicator of radiation or scattering problem',/ 

lOx, ' 1 - Implies radiation problem',/ 
lOx, ' 2 - Implies scattering problem',/ 

lOx, '-2 - Implies scattering including incident wave') 
Variable = 'Radiation or Scattering Problem' 

Call Ichange* Variable, IRS12, 0) 

J = Ir s 12 

if* j. ne. 1 . or. j. ne. 2 . or. j. ne. -2 ) then 

Write**,*) 'Input out of range. Try again' 

Goto 730 
End if 

Write**, 740 ) 

740 Format*' Angle for Computing Radiated or Scattered Field ', 

'(Deg)') 

Variable = 'Theta' 

Call Rchange* Variable, Thrd, 1) 

Var iable = 'Phi' 
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Call Rchanqei Variable. Phrdi 1) 


Write(*. *) 'Angle of Indicent Plane l-Jave of Scattering (Deg)' 
Variable = 'Theta' 

Call Rchange< Variable. Thri. 1) 

Variable = 'Phi ' 

Call Rchanqe( Variable. Phri. 1) 

Else 

INDZI « 0 
WRITE(*. 80) 

80 FORMAT (/IX. 'SPECIFY TYPE OF PATTERN CUT'/iOX. 

'1 GREAT-CIRCLE CUT (VARIABLE THETA) '/10X, 

'2 CONICAL CUT (VARIABLE PHI)'/) 

READ ( *, * ) I CUT 
WRITE (*.90) 

90 FORMAT (/IX. 'FAR-ZONE ANTENNA PATTERN? (Y/N)') 

READ< *, 50) FARZN 
WRITE ( *, 1 10) 

110 FORMAT (/IX, 'VALUE OF CONSTANT ANGLE: ') 

READ ( *, * ) CONANG 
WRITE ( *, 130) 

130 FORMAT(/ IX. 'ANGLE INCREMENT?') 

READ ( *, * ) ANGINC 
WRITE ( *> 120) 

150 FORMAT (/IX, 'DO YOU WANT AN OUTPUT FOR PLOTTING? (Y/N)') 

READ (*,50) PLTOUT 

IF ( FARZN .EG. 'N' OR. FARZN . EQ. 'n ' ) Then 

PARAMETER VALUES SET FOR FAR-ZONE RADIATION 
PATTERN FOR BOTH TYPES OF PATTERN CUT 

Variable = 'Radius for Field Point (meters)' 

Call RChange( Variable, RF» 1) 

I FE=0 
IPFE=0 
FNDFE= i 
PHFE=1. O 
IFA=0 
IPFA=0 
FNDFA=1 
THFA=1 . 0 

INPUT MATRIX DOES NOT EXIST FOR THIS CASE 

I WRZT=0 

PLATES ARE MODELED AS PERFECT CONDUCTORS 

CMM=-1 . 0 

"WIRE RADIUS" SET TO .001 FOR PLATE 

A=. 001 

INFORMATION INPUT FOR PLANE-WAVE SCATTER 

WR I TE ( *, 160) 

160 FORMAT (/IX, 'SPECIFY TYPE OF SCATTER COMPUTATION: '/10X, 

'1 BACKSCATTER '/10X, 

'2 BISTATIC SCATTER '/10X, '3 FORWARD SCATTER') 


C 


32 



o o o o o o o o o o o o o o o o o o 


c 


170 


READ ( *, * ) I TYPE 
IF ( ICUT. EQ. 1 ) THEN 
I SE= I TYPE 
ISA=0 
ELSE 

I SE=0 
ISA=ITYPE 
ENDIF 

IF (PLTOUT. EQ. 'Y,'. OR. PLTOUT. EQ. 'y') THEN 
IPSE=1 
I PSA= 1 
ELSE 

IPSE=0 

IPSA=0 

ENDIF 

- INPUT ANGLE INCREMENTS FOR PLANE-WAVE SCATTER 

FNDSE=ANGINC 

FND5A=ANGINC 

- INPUT CONSTANT ANGLE FOR PLANE-WAVE SCATTER 

PH5E=C0NANG 

THSA=CONANG 

INPUT ANGLE OF INCIDENCE 

WRITE!*, 170) 

FORMAT! /IX, 'WHAT IS THE ANGLE OF INCIDENCE? (THETA, 
READ!*,*) THIN, PHIN 

Else 

PARAMETERS SET FOR PLANE-WAVE SCATTER CASE 

RF = -1 
I SE=0 
IPSE=0 
FNDSE=1 
PHSE= 1 . 0 
THIN=1 . 0 
PHIN=1. 0 
I SA=0 
IPSA=0 
FNDSA=1 
THSA= 1 

- CONSTANT ANGLE VALUE IS SET FOR EITHER PATTERN CUT — 


PHFE=CONANG 

THFA=CONANG 

TYPE OF PATTERN CUT IS SET 

IF! ICUT. EQ. 1 ) THEN 
IFE = 1 
I FA=0 
ELSE 


PHI) ') 
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IFE=0 

IFA=1 

ENDIF 

OUTPUT DESIGNATION ANGLE INCREMENT VALUE — 

IF ( PLTOUT. EQ. ' Y ' . OR. PLTOUT. EQ. 'y'> THEN 
IPFE=1 
IPFA=1 
ELSE 

IPFE=0 

IPFA=0 

ENDIF 

FNDFE=ANGINC 

FNDFA=ANGINC 

End if 
End if 


INPUT FREQUENCY 


175 WRITE**, 180) 

100 FORMAT* IX, 'FREQUENCY IN MEGAHERTZ?') 

READ**,*) FMC 

INPUT WIRE CONDUCTIVITY 8< RADIUS IF WIRES USED 

IF (USEWRS. EQ. 'N ' . OR. USEWRS. EQ. 'n') then 


SUPPLY FICTITIOUS NUMBERS FOR WIRE PARAMETERS - 

210 CMM=-1 . 0 

A=0. 001 

Else 

WRITE**, 190) 

190 FORMAT* IX, 'WIRE CONDUCTIVITY IN MEGAMHOS/METER?') 
READ**,*) CMM 
WRITE**, 200) 

200 FORMAT* IX, 'WIRE RADIUS IN METERS?') 

READ**,*) A 
End if 

215 IF* USEPLT . EQ. 'Y' .OR. USEPLT . EQ. 'y'> Then 

INPUT PLATE INFORMATION 

WRITE**, 220) 

220 FORMAT ( / 1 X , 'HOW MANY PLATES?') 

READ**,*) NPLTS 

WRITE**, 225) 

225 FORMAT* /IX, 'INPUT TYPE OF PLATE TEST MODE: '/10 X, 

'0 FULL-SURFACE PATCH MODE'/lOX, 

'1 FILAMENTARY TEST MODE') 

READ**,*) IFIL 


WRITE**, 230) 
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230 FORMAT (/IX* 'PLATE INFORMATION — V) 

C 

DO 420 NPL=1, NPLTS 

CALL PLATE < NPL, NCNRS, SEGM, IREC, IPN, IGS, PCN, ICHG, 
Zsht) 

430 CONTINUE 

End if 

C 

REUSE OF IMPEDANCE MATRIX 

Disk storage of impedance is not anticipated* therefore: 
IWRZM=0 
IRDZM=0 


IF < USEWRS . EQ. 'Y' OR. USEWRS . EQ. 'yd Then 


INPUT WIRE INFORMATION 

CALL WIRE(NM» NP, NAT, NFPT, NSF1, NSF2, ICHG) 

CALL WRPTS( I, NP. X, Y, Z, ICHG) 

CALL ENDPT (NM, I A, IB, ICHG) 

CALL FPPLT (NFPT, IFM, IABFP, VLG, ZL, ICHG) 

C 

DO 700 1 = 1. NAT 

CALL FPNPL( I, NAS, IABAP, NPLA, VGA, ZLDA, BDSK, ICHG) 
700 CONTINUE 
End if 

c 

C PARAMETER VALUES WRITTEN TO DATA FILE 


C 


810 WRITE (10* 100) NGO, IMPRINT, NRUNS, NWGS, IWR, IWRZT, INT, 

INTD, INWR, IRGM, IFIL, Rf, Indzi 
100 FORMAT (212, 214, 212, 314, 312, F6. 2, 14) 

If< Indzi . ne. 0) Then 

Write(10, *) Fmcl, Fmc2, Dfzi, Dff, Irsl2, Thrd, Phrd, 
Thine, Phinc 


910 

920 

C 


End if 


If ( Indzi . eq. 0) Then 

WRITEdO, 910) IFE, IPFE, 
WRITE( 10, 910) IFA, IPFA, 
WRITE( 10, 910) ISE, IPSE, 
WRITE< 10, 910) ISA, IPSA, 
WRITE* 10, 920) FMC, CMM, 


FNDFE, 

FNDFA, 

FNDSE, 

FNDSA, 

A 


End if 

Format (212, 4F6. 2) 
Format (2F8. 2, F0. 4) 


PHFE 

THFA 

PHSE, 

THSA 


THIN, PHIN 


IF <USEPLT. EQ. 'N'. OR. USEPLT. EQ. 'n') Then 


WRITE( 10, # ) 0 

WRITE< 10, *) 4, 0.2, 1, 3, 0 

DO 1 = 1, 4 

WRITEdO, *) 1.0, 1.0, 1.0 


End do 


Else 

WRITEdO, *) NPLTS 


INTP, 


C 


930 


DO NPL=1, NPLTS 

WRITEdO, 930) NCNRS(NPL) , SEGM(NPL), IREC(NPL), IPN(NPL), 
IGS(NPL) , zsht(Npl) 

Format (13, F6. 3, 313, ' ( ', F6. 2, ', 7 , F6. 2. ' ) ') 

DO NCNR=1, NCNRS(VJPL) 


35 



» 

WRITE < 10, 940) PCN ( 1 , NCNR, NPL ) , PCN ( 2, NCNR. NPL ) , 
PCN (3, NCNR, NPL) 

940 Format (3F7. 2) 

End do 
End do 
End if 


C INPUT FILLER NUMBERS FOR NO-PLATES CASE 

C 

C 

WRITE! 10, *) IWRZM, IRDZM 




950 


960 


- DO NOT WRITE WIRE INFORMATION TO FILE FOR NO-WIRES CASE 

If( Inwr . ne. 0 and. Irgm . ne. 0) Then 

WRITE! 10, 950) NM, NP, NAT, NFPT, NFS1, NFS2 
Format!6I3) 

DO 1 = 1, NP 

WRITE! 10, 960) X ( I ) , Y ( I ) , Z ( I ) 

End do 

Format ( 3FS. 3 ) 


970 

C 

c 


900 

C 

C 

c 

c 

c 


c 


DO 1=1, NM 

WRITE! 10, 950) IA!I),IB!I) 

End do 

IF! NFPT . GT. 0) Then 
DO 1=1, NFPT 

WRITE! 10, 970) IFM ! I ) , IABFP! I ), VLG ! I > , ZL ( I ) 

End do 
Else 

WRITE! 10, 970) 1 , 0, ( 1 . , 0. ) , ( 1 . , 0. ) 

End if 

Format!2I3, ' ! ' , F6. 2, ' , ' , F6. 2. ' ) ( ' , F6. 2, ' , ' , F6. 2, ' ) J ) 


DO 1=1, NAT 

WRITE! 10, 900) NAS! I ) , IABAP ( I ) , NPLA! I ) , VGA ! I ) , ZLDA ! I ) , BDSK! I ) 
End do 

Format!3I3, ' ( ' , F6. 2, ' , ' , F6. 2, ' ) ' , ' ( ' , F6. 2, ' , 1 , F6. 2, ' ) ',F8.3) 

End if , 

END 


Subroutine Rchange ( Variable, Value, New) 


Charac ter-K-1 
Ch arac ter-«-40 
Rea 1 -»-4 
Integer#4 


Change 
Var iab le 
Value 
New 


n = 1 

Do while! Variable 
n = n + 1 
End do 


n = n - 


( n : n + 1 ) 


. ne. 


' ) 


1 
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c 


if< New . eq. 0) then 

Wr i te ( #, 10 ) Var iab 1 e ( 1 : n ) i Value 
10 Formate 1H$) 'The current value of '»A<n>» ' is ' , FQ. 2) 

Write <*, 30) 

30 Format < 1H$» 'Would you like to change this value? ' ) 

Read(^*50) Change 
Else 

Change = ' y ' 

End if 

If< Change . eq. 'y' . or. Change . eq. ' Y ' ) Then 
Wr i t e ( *» 40 ) Var i ab 1 e ( 1 : n ) 

40 Formate 1H$# 'Please enter the new value for ' » A<n>» ' ') 

READ( *) Value 
End if 


50 

100 


<.. 

c 

c 


c. 


Wr ite < *» iOO) 

Format ( A1 ) 

Formate/) 

Return 

End 


Subroutine Cchange ( Variable) Value> New) 


Character*-1 
Charac ter*40 
Comp 1 e x *8 
Integ er*4 


Change 
Var iab le 
Value 
New 


n = 1 

Do whilee Variable (n:n+l) . ne. 

n = n + 1 
End do 
n = n - 1 


' ) 


if ( New eq. 0) then 

Writee*»10) Var iab le ( 1 : n ) > Value 
10 Formate 1H$, 'The current value of ' > A<n># ' is 

' ( F6. 2) F6. 2, ')') 

* 

Wr i te ( #> 30 ) 

.30 Formate 1H$» 'Would you like to change this value? ' ) 

Read<#»50) Change 
Else 

Change = ' y ' 

End if 

I f < Change . eq. 'y' .or. Change . eq. 'Y') Then 
Write <•«•. 40) Var i ab 1 e ( 1 : n ) 

40 Formate 1H$» 'Please enter the new value for ', A<n>» ' ') 

READ ( *> *) Value 
End if 

Wr i te ( 100 ) 

50 FormateAl) 


i 00 Format £ / ) 

Return 

End 

c 

C 

Subroutine Ichange < Variable^ lvalue^ New) 

c 

Character^l Change 

Charac ter *40 Variable 

Integer-*4 New 

Integer*4 lvalue 

c 

n = 1 

Do wh i 1 e ( Variable ( n : n + 1 ) . ne. ' ' .and. n. 1 1 . 40 ) 

n = n + 1 
End do 
n = n - 1 
c 

if( New .eq. 0) then 

Write(^» 10) Variable( 1; n)^ lvalue 
10 Format ( 1 H$i 'The current value of ',A<n>, ' is ', 15) 

c 

Wr i te ( <-t 30) 

30 Format(lH3L 'Would you like to change this value? ' ) 

Read 50) Change 
Else 

Change = ' i\ ' 

End if 

c 

If( Change eq. ' y ' . or. Change eq. ' Y ' ) Then 
Wr i t e ( ■*, 40 ) Variable(lrn) 

40 Format<lHt, 'Please enter the new value for ',ACn>, ' ' ) 

READ ( * ) lvalue 
End if 

Wr i t e < 100) 

50 Format < A1 ) 
i 00 Forma t ( / ) 

R e turn 
End 
c 
c 

SUBROUTINE PLATE ( NPL, NCNRSi SEGM, IREC, IPN, IGS, PCN, ICHG, 

Zsht) 

c 

DIMENSION NCNRS ( 10), SEGM ( 10), IREC( 10), IGS< iO), IPN( 10), PCN<3, 10, 10) 
Comp lex-* 8 Zsh t < 10) 

Charac ter *40 Var iab 1 e 

t: 

WR ITE ( *, 240 ) NPL 

240 FORMAT (IX, 'For PLATE NUMBER ',13) 

c 

Variable = 'Number of Corners' 

CALL ICHaNGE ( Variable, NCNRS (NPL) , ICHG) 

Variable = 'Size of Patch' 

CALL RCHaNGE( Variable, SEGM (NPL) , ICHG) 


WRITE(*, 254) 
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254 FORMAT ( i X* 'IDENTIFY SHAPE OF PLATE -'/ 

10X. 'O POLYGONAL'/ 

10X. '1 RECTANGULAR'/) 

Variable = 'Shape of Plate' 

CALL I CHaNGE ( Variable, IREC(NPL), ICHG) 

r 

C *** POLARIZATION SELECTION IS SET TO 3 FOR ALL CASES *** 

IPN( NPL ) =3 

C 

WRITE!*, 256) 

256 FORMAT! IX, 'WHICH IS THE GENERATING SIDE?') 

Variable = 'The Generating Side' 

CALL ICHaNGE! Variable, IGS!NPL), ICHG) 

C 

Write!*, 500 ) 

500 Format! IX, 'The complex Sheet Impedance must be set to zero ', 

'if the plate contains'/, ' a wire/plate attachment') 
Variable = 'Sheet Impedance' 

Call CChange! Variable, Zsht!Npl), Ichg) 

c 

WRITE!*, 260) 

260 FORMAT ! / 1 X, 'INPUT POSITION OF EACH CORNER (X, Y, Z)'> 

r 

DO 275 NCNR= 1 , NO MRS ( NPL ) 

WRITE!*, 270) NCIMR 
270 FORMAT! IX, 'NO. ',13) 

Variable = 'X coordinate' 

CALL RCHaNGE! Variable, PCN ! 1 , NCNR, NPL ) , Ichg) 

Variable = 'Y coordinate' 

Call Rchange! Variable, PCN!2» NCNR, NPL) , Ichg) 

Variable = 'Z coordinate' 

Call Rchange! Variable, PCN ( 3, NCNR , NPL ) , ICHG) 

275 CONTINUE 
RETURN 
END 


SUBROUTINE WIRE ! NM, NP, NAT, NFPT , NSF1, NSF2, ICHG) 

Character*40 Variable 
Character*l Mcoupl 

Variable = 'the Number of Wire Segments' 

CALL ICHaNGE! Variab 1 e, NM, ICHG) 

Variable = 'Total No. of Points in Wire Structure' 
CALL I CHaNGE (Variable, NP, ICHG) 

Variable = 'Number of Wire-to-Plate attachment pts. ' 
CALL ICHaNGE! Variab le, NAT, ICHG) 

Variable = 'Number of Feed Pts. in Wire Structure' 
CALL ICHaNGE (Variab 1 e, NFPT, ICHG) 

WRITE!*, 470) 

470 FORMAT! 1H*. 'IS MUTUAL COUPLING TO BE COMPUTED?') 
READ!*, 50) MCOUPL 
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50 


FORMAT (A1 ) 

IF ( MCOUPL . EG. 
Variable = '1 
CALL ICHaNGE! 

Variable = '1 
CALL ICHaNGE ( 


' Y " .OR. 
ocation 
Var i ab 1 e 

ocation 
Var i ab 1 e 


MCOUPL . EQ. ' 
of First Feed 
, NSFi. ICHG) 

of Second Feed 
, IMSF2, ICHG) 


NSF.l=0 


NSF2=0 
End if 
RETURN 
END 


y ' ) Then 
Por t " 

Port " 


SUBROUTINE WRPTS! I, IMP, X, Y, Z, ICHG) 
DIMENSION X! 10), Y! 10), Z! 10) 

Character*40 Variable 


WRITE ( *, 520) 

520 FORMAT (/IX, " INPUT LOCATIONS <X, Y, Z) OF WIRE POINTS: 
DO 530 1=1, NP 

WRITE!*-, 525) I 


FORMAT! IX, 'For POINT 

/ 

l 

13) 



Variable = 'X 

coord i 

na t 

e ' 



Call 

Rchange ! 

Var iab 

1 e, 

X! 

i ) i 

ichg ) 

Variable = 'Y 

c oord i 

na t 

e ' 



Call 

Rchange! 

Var iab 

le# 

Y ( 

i ) i 

ichg ) 

Variable = "Z 

c oord i 

na t 

e ' 



CALL 

RCHaMGE ! 

Var iab 

lei 

Z< 

I), 

ICHG) 


530 CONTINUE 
RETURN 
END 


) 


SUBROUTINE ENUPT ( NM, I A, IB, ICHG) 
DIMENSION IA< 10), IB< 10) 


Charac ter*-40 Variable 


WRITE!*, 535) 

535 FORMAT! /IX, ' INPUT ENDPOINT A & ENDPOINT B OF EACH WIRE SEGMENT 
DO 550 J= 1 , NM 


540 

WRITE!*, 540) J 
FORMAT! IX, 'For SEGMENT 

", 13) 



Variable = "Endpoint A" 
CALL ICHaNGE! Variable, 

I A ( J > , 

ICHG) 


Variable = "Endpoint B' 
Coll Ichange! Variable, 

IB! j), 

Ichg ) 


550 CONTINUE 
RETURN 
END 


SUBROUTINE FPPLT ( NFPT, IFM, IABFP, VLG, ZL, ICHG) 
DIMENSION IFM! 10) , 1ABFP!10), VLG!10), ZL(10) 


COMPLEX VLG 
COMPLEX ZL 


Charac ter*40 Variable 
c 

WRITE ( *, 560) 

560 FORMAT (/IX* 'INPUT LOCATIONS OF FEED POINTS: '/> 

DO 650 1=1 . NFPT 
WRITE (*,570) I 

570 FORMAT (IX, 'FEED POINT NO. ',13, ' is in which Wire Segment?') 
Variable = 'Wire Segment' 

CALL I CHaNGE ( Variable, IFM(I), ICHG) 

WRITE(*. 580) 

580 FORMAT (IX, 'AT WHICH END OF SEGMENT? '/ 10X, 

'0 ENDPOINT A'/lOX, 

'1 ENDPOINT B'/) 

Variable = 'Endpoint' 

CALL ICHaNGE ( Variable, IABFP(I), ICHG) 

r 

Variable = 'Generator Complex Voltage at Feed Point' 

Call Cchange( Variable, VLG(I), Ichg) 

c 

Variable = 'Complex Load Impedance (ohms)' 

Call Cchange( Variable, ZL(I), Ichg) 

c 

650 CONTINUE 
RETURN 
END 

c 

c 

C 

SUBROUTINE FPNPL ( NAT , NAS, IABAP, NPLA, VGA, ZLDA, BDSK, ICHG) 
DIMENSION NAS (10), I AB AP ( 10 ) , NPLA ( 1 0 ) , BDSK ( 1 0 ) 

c 

Complex Zlda(l), Vga(l) 

Character*40 Variable 

c 

IF (NAT. EG. 0) Return 

r 

DO 700 1=1, NAT 

Variable = 'No. of wire segments attached to plate' 

CALL ICHaNGE ( Variable, NAS(I), ICHG) 
c 

WR ITE ( *, 660) 

660 FORMAT (/IX, 'ENDPOINT ATTACHED TO PLATE: '/10X, 

'0 POINT A'/lOX, 

'1 POINT B') 

Variable = 'Endpoint' 

CALL ICHaNGE ( Variable, IABAP(I), ICHG) 

c 

Variable = 'Plate Number' 

CALL ICHaNGE ( Variable, NPLA(I), ICHG) 

c 

Variable = 'Attachment Pt Complex Generator Voltage' 

Call Cchange( Variable, VGA(I), Ichg) 

c 

Variable = 'Attachment Point Complex Load Impedance' 

Call Cchange( Variable, ZLDA(I), Ichg) 

c 

Variable = 'Disk Radius (meters)' 

41 



Call Rchange< Mariablei Bdsk(i)< Ichg) 
700 CONTINUE 
RETURN 
END 


A 
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APPENDIX B 
GRAPHICS CODES 
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PROGRAM PLTGTD 
IMPLICIT INTEGER*2(I-N) 

DIMENSION ANG{ 600 ) 

DIMENSION AMP ( 600 ) 

DIMENSION ASTEP ( 5 ) 

DIMENSION XPLT (361,7) 

CHARACTER *1 Z , CHCUT, HOLDIT, NEXTP 
C 

OPEN ( 2, FILE = ' PRN ' ) 

C 

WRITE!*, 70) 

70 FORMAT! IX, 'THE PATTERN IS READY TO BE PLOTTED'//) 

C 

C INPUT THE TYPE OF POLARIZATION COMPONENT DESIRED 

C 

1000 WRITE!*, 100) 

100 FORMAT! IX, 'ENTER POLARIZATION COMPONENT DESIRED :'/ 5X, ’ 1 - E-THETA 
. V5X, '2 - E-PHI ',/5X, '3 - MAJOR AXIS OF POLARIZATION ELLIPSE'/SX 
.,'4 - MINOR AXIS OF POLARIZATION ELLIPSE ' /5X, ' 5 - TOTAL AMPLITUD 
. E 1 ) 

READ!*,*) I POL 
C 

C ASSOCIATE DATA COLUMNS WITH POLARIZATION DESIGNATIONS 

C 

IF(IPOL.EQ.l) KPOL =6 
IF! IPOL.EQ. 2 ) KPOL=7 
IF! I POL . EQ . 3 ) KPOL=3 
IF! I POL . EQ . 4 ) KPOL=4 
IF! IPOL.EQ. 5) KPOL=5 
C 

JCL1=68 

JCL2=566 

JROW1=30 

JROW2=150 

C 

MI NOR=l 
LABEL-1 
NDEC=0 
IOPT=0 
C 

C READ DATA FROM FILE PLOT.DAT TO ARRAY XPLT 

C 

NUMPNT=0 

OPEN (UNIT* 1, FILE* ' PLOT . DAT ' , STATUS® ' OLD ' ) 

DO 150 1*1,361 

READ! 1, 120, END*200 ) XPLT! I , 1 ) , XPLT! I , 2 ) , XPLT! I , 3 ) , XPLT! I , 4 ) ,XPLT( I 
2,5), XPLT (1,6), XPLT (1,7) 

120 FORMAT! 1X,5(3X,F7.2),2(6X,F7.2) ) 

C 

C COUNT POINTS AS READING PROCEEDS 

NUMPNT*NUMPNT+1 
C 

150 CONTINUE 

200 CONTINUE 44 

C 



C DETERMINE TYPE OF PATTERN CUT (GREAT CIRCLE OR CONICAL) 

C CHCUT= 1 Y ' FOR GREAT-CIRCLE CUT; CHCUT='N' FOR CONICAL CUT 

C 

VAL1=XPLT(2,1) 

VAL2=XPLT(1 / 1) 

CHCUT= ' N ' 

IF( ABS ( VAL2-VAL1 ) .GT.0.01) CHCUT='Y' 

C 

C SET MINIMUM AND MAXIMUM VALUES OF ANGLE 

C 

IF(CHCUT.EQ. ' Y' ) THEN 
AFST=XPLT (1,1) 

ALST= ( XPLT( 3,1 ) -XPLT( 2 , 1 ) ) * ( NUMPNT-I ) +XPLT( 1, 1 ) 

NCUT=1 

ENDIF 

C 

I F ( CHCUT . EQ . ' N ' ) THEN 
AFST=XPLT (1,2) 

ALST=XPLT ( NUMPNT, 2 ) 

NCUT=2 

ENDIF 

C 

ANGDI F=XPLT ( 2 # NCUT ) -XPLT ( 1 , NCUT ) 

C 

C TRANSFER DATA FROM XPLT ARRAY TO ANG AND AMP ARRAYS 

C 

ANG( 0 ) =0 . 0 
C 

DO 80 1=1/ NUMPNT 
ANG ( I )=XPLT( I/NCUT) 

IF(ANG(I) .LT.ANG(I-l) ) THEN 
ANG( I ) =ANG( I -1 ) + ANGDI F 
ENDIF 

AMP ( I ) =XPLT ( I / KPOL ) 

IF(AMP( I ) .LT. -50 . ) AMP ( I ) =-50 . 

80 CONTINUE 
C 

C FIND AMPMIN AND AMPMAX VALUES - 

C 

AMPMIN=-10 . 

AMPMAX=10 . 

C 

DO 20 1=1/ NUMPNT 

I F ( AMP ( I ) . GT . AMPMAX ) AMPMAX= AMP ( I ) 

IF(AMP(I ) .LT. AMPMIN) AMPMIN=AMP ( I ) 

20 CONTINUE 
C 

25 ASPECT= ( ALST- AFST ) / ( AMPMAX-AMPMIN ) 

C 

C ROUND OFF AMPMIN & AMPMAX VALUES - 

C 

AMPMIN=INT( AMPMIN) 

AMPMAX* I NT ( AMPMAX ) 

C 

C DETERMINE BEST ANGLE AND AMPLITUDE INCREMENTS TO USE - 


45 



c 


XINC= ( ALST-AFST ) /10 . 
YINC= ( AMPMAX-AMPMIN )/10 . 


ASTEP ( 1 ) =1 . 

ASTEP( 2) =2 . 

ASTEP ( 3 ) =5 . 

ASTEP ( 4) =10. 

ASTEP ( 5 ) =30 . 

C 

DO 30 1=1,5 
XSTEP=ASTEP ( I ) 

IF ( XSTEP . GT . XI NC ) GOTO 40 
30 CONTINUE 
40 XMAJOR=XSTEP 
C 

DO 50 1=1,4 
YSTEP=ASTEP( I ) 

I F ( YSTEP.GT. YINC) GOTO 55 
50 CONTINUE 
55 YMAJOR=YSTEP 

I STEP = I NT ( YSTEP ) 

C 

NSTEP=AMPMIN/ 1 STEP 
IF( AMPMIN . LT . -50 . ) GOTO 58 
AMPMIN=NSTEP* I STEP-I STEP 
GOTO 350 

58 AMPMIN=NSTEP*ISTEP 
C 

350 CALL QSETUP (0,3, -2,0) 

C 

60 CALL QPLOT ( JCLl , JCL2, JROW1 , JROW2, AFST, ALST, AMPMIN, AMPMAX, AFST, AMPM 
1IN,I0PT,1., ASPECT) 

CALL QSMODE ( 6 ) 

C 

YMAX=AMPMAX+0 . 1* ( AMPMAX-AMPMIN ) 

CALL QXAXIS ( AFST, ALST, XMAJOR, MINOR, LABEL, NDEC ) 

CALL Q YAX I S ( AMPM I N , YMAX , YMA JOR , MINOR , LABEL , NDEC ) 

I F ( CHCUT . EQ . ' Y ' ) GOTO 300 

CALL QPTXTA(12, 'PHI, DEGREES ', 2 ) 

GOTO 400 

300 CALL QPTXTA( 14, 'THETA, DEGREES ', 2 ) 

400 CALL QPTXTC ( 4 , ' DB ' , 2 ) 

C 

Jl = l 
C 

CALL QTABL ( J1 , NUMPNT, ANG, AMP ) 

C 

CLOSE ( 1 ) 

CLOSE ( 2 ) 

C 

READ( *, 600 ) HOLDIT 
600 FORMAT (Al) 

C 


WRITE( *, 700 ) 
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700 FORMAT ( //////////////////// ) 

WRITE ( * / 7 5 ) 

75 FORMAT* IX, ' DO YOU WANT TO MAKE ANOTHER PLOT? 
READ( *, 600 ) NEXTP 

IF(NEXTP.EQ. ' Y ' .OR. NEXTP . EQ y ' ) GOTO 1000 

STOP 

END 


Y OR N' 
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PROGRAM PLTMM 
IMPLICIT INTEGER *2 ( I -N) 

DIMENSION ANG( 600 ) 

DIMENSION AMP ( 600 ) 

DIMENSION ASTEP ( 5 ) 

DIMENSION XPLT(361,7) 

CHARACTER *1 Z,HOLDIT,NEXTP 
C 

OPEN (2,FILE=' PRN 1 ) 

WRITE ( * , 70 ) 

70 FORMAT (IX, 'THE PATTERN IS READY TO BE PLOTTED 1 //) 

INPUT THE TYPE OF POLARIZATION COMPONENT DESIRED 
1000 WRITE ( * , 100 ) 

100 FORMAT ( IX, 'ENTER POLARIZATION COMPONENT DESIRED :• /5X, ' 1 - E-THETA 
. '/5X, ' 2 - E-PHI '/) 

READ( *, * ) I POL 

ASSOCIATE DATA COLUMNS WITH POLARIZATION DESIGNATIONS 

IF(IPOL.EQ.l) KPOL=2 
IF( IPOL.EQ.2) KPOL= 4 

JCL1=68 

JCL2=566 

JROW1=30 

JROW2=150 

MINOR=l 

LABEL-1 

NDEC=0 

IOPT=0 

READ DATA FROM FILE PLOT.DAT TO ARRAY XPLT 
NUMPNT=0 

OPEN ( UNIT=1,FILE= ' PLOT . DAT ' , STATUS = ' OLD ' ) 

DO 150 1=1,361 

READ (1,120, END= 200 ) JANG, XPLT (1,2), XPLT (1,3), XPLT (1,4), XPLT (1,5) 
120 FORMAT ( I4,3(1X,F5.1,1X,F6.1) ) 

XPLT (1,1) =REAL ( JANG ) 

COUNT POINTS AS READING PROCEEDS 
NUMPNT=NUMPNT+1 

150 CONTINUE 
200 CONTINUE 

SET MINIMUM AND MAXIMUM VALUES OF ANGLE 
AFST=XPLT ( 1 1) 

ALST= ( XPLT (3,1) -XPLT (2,1)) *REAL ( NUMPNT-1 ) +XPLT (1,1) 

NCUT=1 
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ADIF=XPLT( 2,1)-XPLT(1,1) 

C 

C TRANSFER DATA FROM XPLT ARRAY TO ANG AND AMP ARRAYS 

C 

ANG( 0 ) =0 . 0 
C 

DO 80 1=1, NUMPNT 
ANG(I)=XPLT(I,1) 

IF ( ANG( I ) . LT . ANG( I -1 ) ) THEN 
ANG( I ) =ANG( I -1 ) +ADIF 
ENDIF 

AMP ( I ) =XPLT ( I , KPOL ) 

IF ( AMP ( I ) . LT . -50 . ) AMP (I ) =-50 . 

80 CONTINUE 

FIND AMPMIN AND AMPMAX VALUES - 

AMPMIN=-10 . 

AMPMAX= 1 0 . 

DO 20 1=1, NUMPNT 

I F ( AMP ( I ) . GT . AMPMAX ) AMPMAX= AMP ( I ) 

IF ( AMP ( I ) .LT. AMPMIN) AMPMI N=AMP ( I ) 

20 CONTINUE 

25 ASPECT= ( ALST-AFST )/ ( AMPMAX-AMPMIN ) 

ROUND OFF AMPMIN & AMPMAX VALUES - 

t 

AMPMIN=INT( AMPMIN) 

AMPMAX= I NT ( AMPMAX ) 

DETERMINE BEST ANGLE AND AMPLITUDE INCREMENTS TO USE - 

XINC= ( ALST-AFS f ) /10 . 

YINC=( AMPMAX-AMPMIN )/10. 

C 

ASTEP ( 1 ) =1 . 

ASTEP ( 2 ) =2 . 

ASTEP ( 3 ) =5 . 

ASTEP ( 4 ) =10 . 

ASTEP( 5 ) =30 . 

C 

DO 30 1=1,5 
XSTEP=ASTEP ( I ) 

I F ( XSTEP . GT . X I NC ) GOTO 40 
30 CONTINUE 
40 XMA J OR = X S TEP 
C 

DO 50 1=1,4 
YSTEP=ASTEP ( I ) 

IF( YSTEP.GT. YINC) GOTO 55 
50 CONTINUE 
55 YMAJOR=YSTEP 

I STEP* I NT ( YSTEP ) 
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c 


NSTEP=AMPMIN/ISTEP 
I F ( AMPMIN . LT . -50 . ) GOTO 58 
AMPMIN=NSTEP*ISTEP-ISTEP 
GOTO 350 

58 AMPMIN=NSTEP*ISTEP 

C 

350 CALL QSETUP(0,3,-2,0) 

C 

60 CALL QPLOT ( JCL1 , JCL2 , JROW1, JROW2 , AFST, ALST, AMPMIN, AMPMAX, AFST, AMPM 
1IN, I OPT, 1 . /ASPECT) 

CALL QSMODE ( 6 ) 

C 

YMAX= AMPMAX +0 . 1* ( AMPMAX-AMPMIN ) 

CALL QXAXIS ( AFST, ALST, XMAJOR, MI NOR , LABEL, NDEC ) 

CALL Q YAX I S ( AMPMI N , YMAX, YMA JOR , MI NOR , LABEL , NDEC ) 

300 CALL QPTXTA( 14, 'ANGLE, DEGREES ', 2 ) 

400 CALL QPTXTC( 4, ' DB * , 2 ) 

C 

Jl = l 
C 

CALL QTABL ( J1 , NUMPNT, ANG, AMP ) 

C 

CLOSE ( 1 ) 

CLOSE ( 2 ) 

C 

READ ( *, 600 ) HOLDIT 
600 FORMAT (Al) 

C 

WRITE( *, 700 ) 

700 FORMAT {/////////////////////) 

WRITE(*, 75) 

75 FORMAT(lX, 'DO YOU WANT TO MAKE ANOTHER PLOT? Y OR N ' ) 

READ( *, 600 ) NEXTP 
IF(NEXTP.EQ. ' Y' ) GOTO 1000 
END 
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